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Abstract 

In order to improve the description of the N-d vector analyzing powers A y 
and iT\\ (the so-called A y puzzle) a spin-orbit three-body force has been 
introduced. The L • S term in the NN potential has been modified including 
a two-parameter three-body function depending on the hyperradius p. The 
modification has been performed in channels where the pair of nucleons 
are coupled to spin Sij = 1 and isospin Tij = 1. Calculations have been 
done in the energy region from Ei a b = 648 keV up to E\ a b = 10 MeV. A 
noticeable improvement in the description of the polarization observables has 
been obtained. 

I. INTRODUCTION 

In the standard description of light nuclei as composed of structureless nucleons, three- 
nucleon forces (3NF's) play an important role. The simple assumption of nucleons interacting 
through a pairwise NN potential fails to reproduce the experimental binding energies. This 
suggests that the model has to be extended to include many-body forces in the potential 
energy. First derivations of 3NF's are based on two-pion exchange involving a A excitation. 
The most often used potentials of this kind are the Tucson- Melbourne (TM) M, the Brazil 
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(BR) models, and the Urbana (UR) model which also includes a central phenomeno- 
logical repulsive term with no spin-isospin structure. In practical applications, the chosen 
3NF is adjusted to reproduce the A = 3,4 binding energies. When the calculations are 
extended to describe bound states in p-shell nuclei, a persistent underbinding is observed 
in the mass region with A = 5 — 8 ||. In particular, several excited states are not well 
reproduced indicating that the 3NF could contain a more complicate structure. 

A different problem has been observed in N-d scattering at low energies. The calculated 
vector analyzing powers A y and iTu show an unusual large discrepancy with the experimen- 
tal data HI! . Attempts to improve the description of these observables with the inclusion of 
a two-pion three- nucleon force (27T-3NF) were unsuccessful . This discrepancy has been 
called the A y puzzle since it is, by far, the largest disagreement observed in the theoretical 
description of the three-nucleon system at low energies. Recently, it has been shown that 
the puzzle is not limited to the three-nucleon system but a similar problem appears in the 
calculation of A y in p- 3 He scattering ||. The A y puzzle could be also a signal for different 
forms in the three-body potential |]TTJ,[TT 



In ref. || p-d and n-d scattering have been studied below the deuteron breakup. Using 
the experimental data for cross section and vector and tensor analyzing powers at = 2.5 



and 3.0 MeV from Shimizu et al. ||12|| , it was possible to perform a phase-shift analysis 
(PSA) and to compare the results to the theoretical phases. The conclusion was that small 
differences in the P-wave parameters are responsible for the large disagreement in A y and 
iT\\. It is a general feature that all realistic NN potentials underpredict the splitting in the 
A Pj phases and the magnitude of the €3/2— mixing parameter. When one of the 27r-3NF's is 
included in the Hamiltonian there is not an appreciable reduction of the discrepancy, giving 
in some cases a poorer description. In fact, the inclusion of the 27T-3NF tends to reduce the 
splitting in 4 Pj and to slightly increase e 3 / 2 — . These two opposite effects almost cancel each 
other in the construction of A y and iTu. The operator's form of these particular models of 
3NF's does not include L • S terms which are to a large extent responsible for the splitting 
in the P-waves parameters. 



In the present paper a phenomenological spin-orbit three-nucleon force (SO-3NF) is 
introduced in order to study its effect on the N-d vector analyzing powers A y and iT\\ 
at low energies. The L ■ S term in the NN potential has been modified including a two- 
parameter three-body function depending on the hyperradius p. The two parameters are 
related to the strength and range of the force and they have been fixed with the intention to 
improve the description of A y and iT\\ at E^ & = 3.0 MeV, just below the deuteron breakup. 
Three different sets of parameters have been considered and, accordingly, used to calculate 
scattering observables from Ei a b = 648 keV up to Ei a b = 10 MeV. The SO-3NF has been 
introduced in channels where the pair of nucleons are coupled to spin SV,- = 1 and 

isospin Tij = 1. At the level of the two-nucleon (2N) system this channel is related to 
scattering in odd waves. In N-d scattering the 4 Pj parameters and hence A y and iTu are 
very sensitive to the force in this particular channel. On the other hand its effect on the 
binding energy of the three-nucleon system is very small. This is the opposite behavior from 
that produced by the 27T-3NF which gives the main contribution in the J = l/2 + state and 
has small influence on the vector observables. Therefore these two different classes of 3NF's 
are to some extent complementary and can be studied separately. 

The calculations of the N-d scattering observables presented here in the following 
at different energies have been done using the Pair Correlated Hyperspherical Harmonic 
(PHH) JL3| basis. In this method the wave function of the system is expanded in terms 
of correlated basis elements and the description of the system proceeds via a variational 
principle. Bound states are obtained using the Rayleigh-Ritz variational principle whereas 
scattering states are obtained using the generalized Kohn variational principle. This tech- 
nique has been extensively discussed in refs. [f7t |I^|15| for energies below the deuteron breakup 



threshold and, very recently, also applied to energies above the breakup threshold p||T7|l . 

The paper is organized as follows. In Sec. II the two parameter spin-orbit three-body 
force is introduced. In Sec. Ill the polarization observables as well as the binding energy 
of 3 He are studied for specific values of these parameters. The conclusions and perspectives 
are given in the last section. 



II. A TWO-PARAMETER SPIN ORBIT THREE-BODY FORCE 



Disregarding for the moment the presence of three-nucleon forces, the potential energy 
operator of the three-nucleon system is 

^v = £vw(U), (i) 

i<j 

where V-zn is the NN interaction that, in general, is constructed by fitting the 2N scattering 
data and the deuteron properties. Recently, several potential models have been determined 
including explicitly charge dependence which describe the 2N data with a x 2 P er datum « 1. 
Here we will refer to the Argonne AV18 interaction, which is one of these new generation 



potentials [nj|. The nuclear part of the AV18 potential consists of a sum over 18 different 
terms. The first 14 terms are charge independent whereas the four additional operators 
introduce charge symmetry breaking. Each of the first 14 terms includes a projector Pst^j) 
onto the spin-isospin states 5", T of particles (i, j) multiplied by one of the following operators, 
O p = 1, S 12 , L • S, L 2 , (L • S) 2 . The strength of each term is given by a scalar function v v ST {r ij) 
depending on the relative distance between particles For example the L • S interaction 

between particles is defined in the two channels with isospin Tij = 0, 1 and spin = 1 
and is given by the functions i4o( r y) and t4i( r ij)j respectively. 

In the three-nucleon system odd-parity scattering states are directly related to the po- 
tential in channels where particles are coupled to = 1, = 1. In particular here 
we are interested in the L ■ S term and its relation to the vector observables A y and %T\\. 
Limiting the discussion to the L • S interaction in this particular channel, the corresponding 
term in the three-nucleon potential energy is 

V& Y,^(r, J )L, r S, J l><:( U ) . (2) 

i<j 

At this point we can conjecture that the presence of the third nucleon (particle k) could 
modify the interaction, introducing in the function v[\ a dependence on the distances between 
the third nucleon k and particles Accordingly, the above interaction transforms to 
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V 3N = E o[ u, n( r «*) L « ' S 'J + ' Si^n(r ijfc )]P n (zj) , (3) 

i<j 

where is a scalar function of the three interparticle distances r^-, r^, r^. The symmetric 
form has been introduced since, in general, the L • S operator does not commute with 
an operator depending on r^. Different forms are possible for the three-body interaction 
w ii( r ijk) provided that lUiifajfc) ~~ ► u u( r ij) when r^r^ — > oo. A simple two-parameter 
form that we are going to analyze here is 

w l 1 s 1 (r ljk ) = v[\(r tJ ) + W e-^ , (4) 

where the hyperradius p is 

P 2 = ^(r 2 i2 + r 2 23 +r 2 31 ) (5) 

and Wq and a are parameters characterizing the strength and range of the three-body 
term. When the dependence in the scalar function is limited to and p, the operators 
w ii( r ijiP) anc ^ ' commute. Accordingly, the spin-orbit force becomes 

V& = $>i.i(r«)L - • SyPiiCy) + Woe-^^H- ■ SyPnCij) • (6) 

Replacing this term in eq.(^) and including now also the 27T-3NF, the final form for the 
three-nucleon potential energy to be used in the present work is 

VW = £VW(i,j)+ E K N (hj,k)+ £ (7) 

i<j i<j<k i<j<k 

The Wjjv term is the phenomenological spin-orbit force defined in the second term of eq. (H) ; 
for the term the discussion will be limited to the Urbana force. 

The choice of the two-parameter exponential form in the definition of u^i (see eq.(|])) 
is arbitrary and was selected in order to make a phenomenological representation of a 3NF 
containing a spin-orbit interaction. The hyperradial dependence is the simplest scalar func- 
tion depending on the three-interparticle distances which has the property of commuting 
with the spin-orbit operator. These choices, driven by simplicity, have been made in order 
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to focus on the spin-orbit operator and its relation to A y and %T\\ in p-d scattering. The 
numerical values for the constants Wo, a are discussed in the next section. 

In principle, the argument used to introduce the scalar function u4i( r ijfc) as a modifi- 
cation of the function t>n( r ij) due to the presence of particle k, could be extended to the 
other functions Vg T (rij) of the NN potential. This leads to a three-body force with an 
hybrid form in which nucleons interact with the same operator's structure of the NN 
potential but with scalar functions depending on the interparticle distances of the three 
nucleons k). The original function v p ST (rij) must be recovered when nucleon k is at oo. 
In the present work we are limiting the argument to one specific term and trying to relate 
its parametrization to the two vector observables A y and iT n . 



III. N-D SCATTERING CALCULATIONS WITH THE SO 3NF 

The inclusion of the SO-3NF has been done with the hope of improving the description 
of A y and iTu without destroying the agreement already observed for the cross section and 
tensor observables. At energies below the deuteron breakup, measurements for tensor and 
vector observables exist for p-d scattering at Ei^ = 648 keV JlO and at E^ = 2.5 MeV and 



3.0 MeV pHH. The 648 keV data set is at the lowest energy at which these measurements 



have been made. The other two sets of measurements lie just below the deuteron breakup. 
Theoretical calculations using the AV18 potential underpredict A y and iTu of about 30%. 
Comparisons to phase-shift and mixing parameters extracted from PSA performed at these 
three energies |||ll],|l9|] show the aforementioned insufficient splitting in the 4 Pj phase shifts 
as well as an underpredicted €3/2—- Calculations using the AV18+UR potential essentially 
do not change these findings. The almost constant underprediction in the vector polarization 
observables (in percentage) and the fact that the inclusion of the 27i~-3NF's does not increase 
the splitting in A Pj is a motivation for considering new additional forms for the three-body 
potential. The selection of the L ■ S operator is a natural choice since when applied to the 
4 Pj state it acts with opposite sign in the states J = l/2~ and J = 5/2~, tending to increase 
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the splitting. 

Three different choices of the exponent a in the hyperradial spin-orbit interaction defined 
in eq.(f|) have been selected with the intention of constructing forces with different ranges. 
The strength Wq has been adjusted in each case in an attempt to improve the description of 
the vector observables. The analysis has been performed at = 3.0 MeV. The selected 
ranges are a = 0.7, 1.2, 1.5 fm -1 , so as to simulate a long, medium and short range force. 
The corresponding values for the depth are Wq = —1, —10, —20 MeV. The calculations have 
been performed using the nuclear part of AV18 plus the Coulomb interaction. The 27T-3NF 
has been disregarded at the present stage since its contribution to the description of the 
vector observables is small. 

The results for the proton and neutron analyzing powers A y and the deuteron analyzing 
power iT\i are given in Fig.l together with the experimental data of ref. fll2 , |20|| . The four 



curves correspond to the AV18 potential and the three different choices for the parameters 
(a, Wq). The dotted line is the AV18 prediction and shows the expected discrepancy. The 
solid line corresponds to the AV18 plus the long range force (AV18+LS1), the long-dashed 
line to the AV18 plus the medium range force (AV18+LS2) and the dotted-dashed line to the 
AV18 plus the short range force (AV18+LS3). The inclusion of the spin-orbit force improves 
the description of the vector observables, although there is a slightly different sensitivity in 
A y and iT n . The AV18+LS1 curve is slightly above the data, especially for iT n . The 
AV18+LS2 curve is slightly below (above) the data in A y (iTn). The AV18+LS3 curve is 
slightly below the data especially for A y . In the bottom panel of Fig.l the n-d analyzing 
power has been calculated using the same potential models as before. Again, there is an 
improvement in the description of A y equivalent to that one obtained in the p-d case. In 
Fig. 2 the tensor analyzing powers T 2 o,T 2 i,T 2 2 are shown at the same energy and compared 
to the data of ref. (T^|. The inclusion of the SO-3NF has no appreciable effect and the 
four curves are practically on top of each other. These observables are not very sensitive to 
the splitting in 4 P/-waves. They are sensitive to scattering in D-waves and higher partial 
waves, which are only weakly distorted by the SO-3NF. 



Before extending the calculations to other energies the analysis of the binding energy of 
3 He deserves some attention. We expect that the inclusion of the W^ S N term will produce 
only a small distortion in the bound state due to the low occupation probability of channels 
with Sij = 1 and = 1. This is corroborated by the calculations shown in table I. 
The binding energy, the kinetic energy and the S'-, D-, and P-wave probabilities are given 
for the different potential models. The AV18 and the AV18+UR interactions have been 
considered. The 2-7T-3NF of Urbana has been now taken into account since it produces large 
effects in the bound state. Calculations have been done for these two potentials with and 
without the inclusion of the three different choices for the SO-3NF. Its inclusion produces 
a small repulsion which in any case is not greater than 50 keV, with all the other mean 
values modified very little. The force with the longest range produces very tiny effects 
due to its very small strength. The other two forces produce slightly greater, and similar, 
modifications. We can conclude that the structure of the three-body bound state remains 
essentially unaffected by the spin-orbit 3NF with the ranges and strengths considered. 

The analysis of the bound state is important since changes in A y and iTn of the size 
shown in Fig.l could, in principle, be obtained with modified forms of the 27r-3NF's. But, 
in general, these modifications are not anymore compatible with a correct description of the 
bound state. This is not the case for the spin-orbit 3NF we are considering. Therefore, we 
can extend the calculations to other energies based in the fact that the new model makes a 
selective and appreciable effect only in the vector observables. Calculations have been done 
at Ei a b = 648 keV and 2.5 MeV. The results for A y and iT n are shown in Fig. 3 and compared 



to the experimental data of refs. [jTOjJT^I . Again a remarkable improvement in the description 
of both observables is obtained. At the lowest energy the AV18+LS1 model seems to be 
more effective because of its long range whereas for the other two ranges the centrifugal 
barrier plays some role. The results at Ei a i = 2.5 MeV have the same characteristics as 
those at 3.0 MeV. 

In table 2 we compare results for the P-waves parameters at P^ = 2.5 and 3.0 MeV. 
Again, different cases using the AVI 8 potential with and without the inclusion of the spin- 
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orbit 3NF have been considered. In the last column the parameters corresponding to the PSA 
from ref. |J are given for the sake of comparison. In the three cases the SO-3NF increases 
the splitting of the 4 Pj phases. The phases 4 -Pi/2 and 4 -Ps/2 and the mixing parameter 
63/2-, which play a major role in the description of the vector observables, are now in much 
better agreement with the PSA values. It was not obvious from the beginning that with two 
parameters (a, Wo), it would be possible to increase the difference AP = 4 P 5 / 2 — 4 -Pi/2 in the 
measure suggested by the PSA, together with a change in e 3 / 2 _ in the expected direction and 
magnitude. At the three energies the change in €3/2- was observed to happen in the correct 
direction. Not all parameters are closer to the PSA values after including the SO-3NF. For 
example 4 P 3 /2 has slightly changed in the opposite direction to that suggested by the PSA. 
But the final result in the description of the observables was found always to produce an 
improvement. 

In order to complete the study of the SO-3NF the extension to p-d calculations above 
the breakup channel is now considered. This extension is not straightforward since when the 
breakup channel is open a correct description of the three outgoing nucleons has to be done. 
In particular, in the p-d case, the Coulomb interaction introduces difficulties that have been, 
and are at present, subject of intense investigation. Recently the PHH technique has been 
extended to describe p-d elastic scattering above the deuteron breakup [|l6i|i7| . The method 



is based on the use of the Kohn variational principle in its complex form JT5[ and it provides 
an accurate description of the polarization observables. In ref. [[H] differential cross sections 
and vector and tensor analyzing powers at = 5 and 10 MeV have been calculated using 
the AV18 interaction. In the present work, calculations have been done at the same two 
energies using the AV18 interaction with and without the SO-3NF. The results for A y and 
iTxi are given in Fig. 4 together with the experimental data from ref. |HJ (Ei^ = 5 MeV) 
and ref. |S] (Ei a h = 10 MeV). Also at these energies we observe the same trend as before, 
both observables are now better described. The sensitivity to the different ranges is slightly 
different with the medium range model starting to be more effective. In fact, at 5 MeV the 
long and medium range curves overlap and at 10 MeV the splitting in the curves is now 
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reversed in that the upper curve corresponds to the medium range force. 

Finally, let us extend the analysis of the effects of the S0-3NF to other observables 
at different energies. In Fig. 5 the effects on the elastic cross section are shown in the 
energy interval from Ei a b = 648 keV up to Ei a b = 10 MeV. At each energy, the four curves 
corresponding to the different potential models are on top of each other. There is not 
enough sensitivity in this observable to changes in the phase-shift and mixing parameters of 
the magnitude introduced by the SO-3NF used. The tensor analyzing powers at Eid = 5 
MeV and 10 MeV are shown in Fig. 6. The effects on these observables are very small. The 
three curves corresponding to the AV18 plus the SO-3NF are practically superposed with 
slight differencies with the AV18 curve (dotted line). In particular, the second minimum of 
T 20 is slightly improved as well as the second maximum of T 2 i, whereas T 22 is essentially 
unchanged. 

In order to give a quantitative measure of the improvement introduced by the SO-3NF in 
the description of the data, a x 2 (P er datum) analysis is displayed in table III. The analysis 
includes the cross section, the two vector and the three tensor analyzing powers at E^ = 648 
keV, 3 MeV, 5 MeV and 10 MeV. The data set has been taken from ref. |njjggl|-[23| . The 



X 2 per datum has been calculated using the theoretical predictions of the AV18 potential 
model with and without the inclusion of the SO-3NF. The cases corresponding to the three 
different choices of the strength and range parameters have been considered. From the table 
is clear the selective effect of the spin-orbit force on the vector observables. There is a 
dramatic improvement in terms of x 2 i n these observables, whos value is reduced by more 
than one order of magnitude in several cases. At each energy, the AV18+LS1 and AV18+LS2 
potential models give the best description. At the first three energies both potential models 
reproduce the data reasonably well and with similar quality. At 10 MeV the obtained x 2 
per datum for A y is very high in all cases. This is a consequence of the extremely small 
error bars in the data set at this particular energy. However, the x 2 has been improved by 
a factor of « 25 going from AV18 to AV18+LS2 and a further reduction could be obtained 
with a fine tune of the strength and range parameters Wq, a. 
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Looking at the differential cross section, in all cases the x 2 per datum is a large number 
and it changes very little when the SO-3NF is included. There is a sensitivity to the 3 He 
binding energy which is not well reproduced unless a 27T-3NF is considered. For example, 
the value x 2 — 30.3 at 3 MeV obtained with the AV18 potential reduces to x 2 — 4.0 with the 
AV18+UR potential. In the case of the tensor observables the changes in terms of x 2 when 
the SO-3NF is considered are moderated. There is a slight improvement in the description 
of T 20 and T 2 i, whereas the reverse situation is seen in T 22 . 

IV. CONCLUSIONS 

Elastic N-d scattering has been studied in the energy range from E^ = 648 keV to 
Eiab = 10 MeV using a potential model which includes a spin-orbit three- nucleon interaction. 
This three-body potential was introduced as a "distortion" of the function v[\(rij) of the 
NN potential. This function gives the magnitude of the L • S interaction in the channels 
where particles are coupled to spin Sij = 1 and isospin = 1 and was converted to a 
function ti4i( r yfc) depending on the three-interparticle distances r^-, r^, r^. The condition 
on u^i was that the interaction v l {[ is recovered when the third particle is far from the other 
two. A phenomenological hyperradial exponential form depending on two parameters was 
used which fixes the range and strength of the three-body part of the interaction. 

The choice of the L • S operator acting on the — 1, = 1 spin-isospin channels was 
based on the large effects it has in the description of the two vector observables A y and iT u 
in N-d scattering. The origin of this discrepancy lies in the too low splitting predicted by all 
realistic potential models for the A Pj phase shifts. There is a close relation, due to the Pauli 
blocking, between scattering in 4 Pj- waves and the interaction in this channel. Moreover the 
L • S operator is attractive (repulsive) in the J = l/2~ state (J = 5/2~ state), therefore 
increasing the splitting of these phases. 

The study of the SO-3NF was done phenomenologically by fixing three different values for 
its range parameter a and selecting the strength parameter W to provide a better description 
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of the vector observables. The energy at Ei ab = 3.0 MeV was chosen for this analysis and 
successively the three sets of parameters were used to describe the same observables at 
different energies. The difference between the curves obtained for the description of A y and 
iT\\ after including the spin-orbit force could be further reduced with a fine tuning of Wq. 
An important check was made resulting in the observation that the structure of the bound 
state and other observables, for example the tensor analyzing powers are not appreciable 
disturbed by this interaction. Moreover the n-d and p-d A y are described equally well when 
the spin-orbit 3NF is included. This means that no additional charge-symmetry breaking 
effects are needed. 

The use of this phenomenological SO-3NF improves the description of the vector ob- 
servables in all the studied cases, from Ei ab = 648 keV to E iab = 10 MeV. The sensitivity 
of the observables to the range parameters is slightly different at the energies studied and 
there is not a definitive preference for one of them. Perhaps the set with the medium range 
parameter (a = 1.2 fm" 1 , W = —10 MeV) gives on average the best description. A bet- 
ter evaluation about the simple hyperradial dependence introduced here because of its nice 
property of commuting with the L • S operator, can be obtained only after extending the 
calculations to higher energies. The main conclusion of the present work is the identification 
of the L • S force in the Sij — 1, T i} ■, — 1 spin-isospin channel as one which can resolve the 
A y puzzle. In addition, a simple model has been proposed to repair the discrepancy. 

Further investigations, which are in progress, consist in the inclusion of the 27T-3NF in 
the description of the scattering observables other than in the bound state, the extension of 
the calculations to higher energies and the study of the force in the four-body reaction p- 
3 He. Whereas in the first case we can expect at most a small change of the two parameters 
(a, Wo)) the two other studies will give more insight about the force. In particular, by 
studying the p- 3 He reaction the extension to heavier systems can be tested. 

Finally, the possibility of introducing a dependence on the three interparticle distances 
in other NN functions Vg T (rij), as discussed at the end of Sec. II, deserves some attention. 
Disregarding other types of 3NF's, the modified potential energy between three nucleons 
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would be 

Vsn = E E E \ Kt(^)^' + 0> P ST(r ijk )]PsT(ij) ■ (8) 

i<j S,T V Z 

The original functions v p ST are recovered when nucleon k is far from the pair In 
addition, the symmetric form can be avoided if the three-body radial dependence is limited 
to the hyperradius, i.e. w p ST {r\j k ) = w p ST {r-ij, p). The next step is the parametrization of 
the functions w p ST with a number of parameters to be determined by a fit procedure. The 
maximum number of free parameters and the selection of the channels where the distorted 
functions are introduced are related to the observables included in the fit. For example, 
it would be impossible to reproduce simultaneously the binding energy of 3 He and the 
vector analyzing powers in p-d scattering by modifying the potential only in channels with 
Sij = l,Tij = 1 and neglecting the 27T-3NF. Conversely, this would be possible by extending 
the modification to other terms in channels with = 1,7^- = or = 0,7^ = 1. 
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Table Captions 

Table 1. The binding energy of 3 He is shown together with the kinetic energy and the 
S'-, P- and .D-wave probabilities. The AV18 and the AV18+UR potentials are considered 
with and without the spin-orbit 3NF for three choices for the parameters (see text). 

Table 2. The P- waves phase shift and mixing parameters calculated at two different 
energies. The AV18 potential and three different sets for the parameters in the spin-orbit 
3NF (see text) have been used. The results from the PSA of ref. || are given in column 6 
for the sake of comparison. 

Table 3. \ 2 P er datum at four different energies calculated using the AV18 potential 
model with and without the three different parametrization of the SO-3NF. The data set 
has been taken from refs. [p!0|, p!2] , pl| -p3l . The number in parenthesis is the number of data 
points. 
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Figure Captions 

Fig.l. The p-d and n-d analyzing power A y and the deuteron analyzing power zT n 
are shown at Eiab = 3.0 MeV. The different curves correspond to the following potential 
models: AV18(dotted line), AV18+LS1 (solid line), AV18+LS2 (dashed line), AV18+LS3 
(dotted-dashed line). The experimental data are from ref. |TIJ (p-d) and ref. |ZD| (n-d). 

Fig. 2. The tensor analyzing powers T20, T21, T22 at E\ a b = 3.0 MeV for the same potential 



model as Fig.l. The experimental data are from ref. [|12 



Fig. 3. A y and %T\\ at two different lab energies. The four curves correspond to the same 
potentials as in the previous figures. Experimental data are from ref. |10| at 648 keV and 
from ref. at 2.5 MeV. 

Fig. 4. As in Fig. 3 at two different energies above the deuteron breakup threshold. Ex- 



perimental data are from ref. |2T] at 5 MeV and from ref. [^] at 10 MeV. 

Fig. 5. The differential cross section at five different lab energies. At each energy the 
four curves corresponding to the AV18 potential with and without the SO-3NF overlap. 
Experimental points are from refs. fT2||2~T| , |2"B[ . 

Fig. 6. The tensor analyzing powers at Ei a b = 5 MeV and 10 MeV. At each energy there 
are four curves corresponding to AV18 (dotted line), AV18+LS1 (solid line), AV18+LS2 
(dashed line) and AV18+LS3 (dotted-dashed line). 
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TABLES 



Potential 


B(MeV) 


T(MeV) 


Ps'(%) 


Pp{%) 


Pd{%) 


AV18 


-6.942 


45.72 


1.542 


0.065 


8.481 


AV18+LS1 


-6.929 


45.66 


1.546 


0.064 


8.456 


AV18+LS2 


-6.905 


45.52 


1.555 


0.063 


8.423 


AV18+LS3 


-6.905 


45.52 


1.556 


0.063 


8.427 


AV18+UR 


-7.768 


50.25 


1.251 


0.132 


9.262 


AV18+UR+LS1 


-7.751 


50.17 


1.255 


0.131 


9.233 


AV18+UR+LS2 


-7.718 


49.99 


1.264 


0.129 


9.190 


AV18+UR+LS3 


-7.718 


49.97 


1.265 


0.129 


9.194 



TABLE I. 
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Ei u — 3 MeV 






AV18 


AV18+LS1 


AV18+LS2 


AV18+LS3 


PSA 


1/2 


-7.36 


-7.40 


-7.38 


-7.37 


-7.41 


* 3/2 


-7.12 


-7.11 


-7.12 


-7.12 


-7.18 


4 p 

1/2 


22.11 


21.67 


21.72 


21.87 


21.77 


*3/2 


24.23 


24.10 


24.14 


24.18 


24.30 


4 P K /o 
M5/2 


24.00 


24.27 


24.25 


24.17 


24.26 


e l/2— 


5.72 


5.62 


5.62 


5.66 


5.70 


o/z — 


-2.22 


-2.37 


-2.33 


-2.29 


-2.46 


K, ,—25 MeV 






AV18 


AV18+LS1 


AV18+LS2 


AV18+LS3 


PSA 


1/2 


-6.85 


-6.88 


-6.87 


-6.86 


-7.11 


2 p 
*3/2 


-6.70 


-6.67 


-6.68 


-6.68 


-6.90 


4 Pl /9 

1/2 


20.10 


19.70 


19.75 


19.87 


19.89 


^3/2 


22.45 


22.27 


22.30 


22.34 


22.49 


^5/2 


21.84 


22.08 


22.07 


22.00 


22.14 


ei/2- 


4.95 


4.87 


4.87 


4.90 


4.77 


e 3/2- 


-1.86 


-2.01 


-1.98 


-1.94 


-2.06 



TABLE II. 
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Einh = 648 keV 


AV18 


AV18+LS1 


AV18+LS2 


AV18+LS3 


a (19) 


38.5 


39.5 


41.0 


41.5 


A, (3) 


22.7 


0.4 


2.6 


8 




1 8 3 

±0.0 


3 4 


1 4 


4 3 


Ton ('24") 


2 7 


2 fi 


2 fi 


2 7 

— ■ • 1 


T 2 i (22) 


3.9 


3.9 


3.9 


3.9 


T 22 (20) 


2.1 


2.1 


2.1 


2.1 


Ei nh = 3 MeV 


AV18 


AV18+LS1 


AV18+LS2 


AV18+LS3 


o- (37) 


30.3 


31.2 


32.9 


32.8 


A, (38) 


236.7 


3.3 


5.8 


46.6 




1 4Q 1 






1 4 R 


Ton f51 


4.4 


2 7 

— ■ • 1 


3 


3 2 


T 2 i (51) 


6.8 


5.0 


4.4 


4.5 


T 22 (51) 


12.7 


14.8 


14.8 


14.3 


E, nb = 5 MeV 


AV18 


AV18+LS1 


AV18+LS2 


AV18+LS3 


a (41) 


29.3 


30.2 


30.2 


30.5 


A, (43) 


229.8 


6.9 


8.5 


38.5 




1 2D 1 


2 7 


^ 2 


5 "\ 
0.0 


Ton ffi'H 


12.2 


8 8 


8 1 


9 1 


T 2 i (61) 


11.2 


9.3 


8.6 


9.2 


T 22 (61) 


7.9 


9.6 


9.7 


9.2 


E, nh = 10 MeV 

J —'L(LD v 


AV18 


AV18+LS1 


AV18+LS2 


AV18+LS3 


o- (42) 


12.3 


11.4 


10.4 


11.1 


A, (48) 


840.9 


80.1 


34.9 


128.3 


iT n (27) 


37.3 


4.2 


3.1 


4.9 


T 20 (27) 


12.8 


8.4 


6.7 


7.9 



20 



T 2 i (22) 


6.3 


3.9 


3.5 


4.0 


T 22 (27) 


25.4 


29.8 


32.0 


29.4 



TABLE III. 
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45 90 135 180 




45 90 135 180 
cm. angle 



45 90 135 180 



cm. angle 




cm. angle cm. angle 




cm. angle cm. angle 



da/da 

[mb/sr] 

300 





da/dQ. 

[mb/sr] 

300 





do/dQ 

[mb/sr] 

300 








da/dQ. 

[mb/sr] 

300 






doldQ 

[mb/sr] 



300 








45 90 135 
cm. angle 



180 



45 90 135 180 
cm. angle 



